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Abstract 

Recently, the possibility was pointed out by one of the present authors and his collaborators 
that an effective naked singularity referred to as "a visible border of spacetime" is generated by 
high-energy particle collision in the context of large extra dimensions or TeV-scalc gravity. In 
this paper, we investigate the particle creation by a naked singularity in general dimensions, while 
adopting a model in which a marginally naked singularity forms in the collapse of a homothetic 
lightlike pressureless fluid. We find that the spectrum deviates from that of Hawking radiation 
due to scattering near the singularity but can be recast in quasi-thermal form. The temperature 
is always higher than that of Hawking radiation of a same-mass black hole, and can be arbitrarily 
high depending on a parameter in the model. This implies that, in principle, the naked singularity 
may be distinguished from a black hole in collider experiments. 
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1 Introduction 



The higher-dimensional scenarios with large pQ and warped [2] extra dimensions were proposed to 
resolve or reformulate the hierarchy between the gravitational and electroweak interactions. In these 
scenarios the d-dimensional {d > 4) Planck mass Mp is related to that in 4-dimension Mp(4)(~ 
10^9 GeV) by Mj,^^^ ~ L'^-^Mp'^, where L is the size of an extra dimension. Thus, if L is large 
enough, the fundamental Planck mass Mp can be as low as a few TeV. For instance, such a small 
Planck mass is realized if L 0.1 cm when d = 6, L ~ 10~^ cm when d = 7. If the Standard Model 
particles but the gravitons (and possibly other unobserved particles) are confined to our 3-brane these 
scenarios are consistent with all current observations. 

One of the most striking predictions of such scenarios is the productions of a large number of mini 
black holes in high-energy particle collisions [3]. A simplified picture of the black hole production 
is put in the following way. The size of a black hole is characterized by Schwarzschild radius rn, 
scaling with its mass as vh oc M^jj^ ^\ If colliding particles (partons) have a center-of-mass energy 
Mbh above a threshold energy of order Mp and an impact parameter less than the Schwarzschild 
radius, a black hole of mass MpH is produced. In other words, the total cross section of black hole 
production is given by apH — ^^h- The black holes so produced will decay thermally via the Hawking 
radiation [Ij and be detected in terrestrial collider experiments such as CERN Large Hadron Collider 
and in ultrahigh energy cosmic rays. Such possibilities have been extensively studied and known to 
give rise to rich phenomenology {e.g., see [5] for a recent review). 

Quite recently, one of the present authors and his collaborators pointed out another interesting 
possibility [6] of the high-energy particle collisions in the TeV-scale gravity. They argued that effective 
naked singularities called the visible borders of spacetime would be generated by high-energy particle 
collisions. A border of spacetime, proposed in [7] originally, is defined as a domain of spacetime where 
the curvature becomes trans-Planckian and acts as an border (or boundary) of classical sub-Planckian 
regimes. A simplified picture of the generation of a visible border is as follow. Suppose colliding 
particles have a center-of-mass energy above Mp, just like the black hole production mentioned above. 
Then, suppose that the impact parameter is small enough to make the energy density of the colliding 
region become trans-Planckian but the impact parameter is larger than the Schwarzschild radius. In 
such a case, if any, the curvature around the colliding region becomes trans-Planckian through the 
Einstein equation, but the horizon will not forir0. Thus, the trans-Planckian domain of spacetime, 

^It was argued in [8], replying to the claim raised in [3, that the curvature singularity at the collision event described 
by the Aichelburg-Sexl waves can be avoided if one takes into account the finite extension of colliding particles as a 
quantum effect. On the other hand, the argument of visible-border production in [B] does not depend on how to describe 
the collision but on a dimensional estimation of the energy density and curvature around the collision event. Furthermore, 
the argument takes into account the possible quantum uncertainty of particles' position. Thus, at this point, we believe 
that the visible-border production would not be avoided even if one takes into account the semiclassical effects such as 
the uncertainty. Of course, however, it would be important to verify further the reality of the visible-border production 
from various semiclassical points of view. This paper, in which another semiclassical effect — the particle creation — 
is investigated, can be regarded as a part of such verifications. The authors thank an anonymous referee for making us 
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i.e., the border of spacetime, is visible or naked to outer observers. Paper [6] showed with a simple 
argument that such phenomena can occur in collider experiments, that is regarded as an effective 
violation of Cosmic Censorship Hypothesis [10]. In such a visible border production, in contrast to 
the black hole production, the trans-Planckian regime is exposed to observers, and therefore an arena 
of quantum gravity could be provided. 

It would be fare to say that the production of visible borders (or naked singularities) as well as the 
black holes in collider experiments is just a possibility: we do not have any rigorous quantum theory 
of gravity to predict phenomena near the Planckian regime. Nevertheless, it would be important to 
predict possible phenomena with using classical/semiclassical tools available. The study of semiclas- 
sical effects during naked singularity formation (in 4 dimensions) has a relatively long history, which 
goes back to the seminal works by Ford and Parker [11] and by Hiscock et al. [12] . (An incomplete list 
of studies in this direction is [131 [H]. See [15] for a review.) Typically, when a (strong) global naked 
singularity forms, which violates even the weak version of the Cosmic Censorship Hypothesis [TO], the 
power of particle creation diverges at the Cauchy horizon (if estimated without backreactions). 

Due to the universality of particle creation by naked singularities, one can expect that similar 
phenomena will occur in higher dimensions, which is addressed in this paper. However, we should note 
several points. Firstly, the generation of visible borders or naked singularities in large extra dimension 
scenarios would be a highly asymmetric phenomenon: the gravity propagates in all dimensions; the 
Standard Model particles are confined to the 3-brane; the non-zero impact parameter is essential. 
It seems difficult to model the visible-border formation by a known solution to Einstein equations. 
Therefore, in this paper, as a first step, we adopt the Vaidya solution describing the spherically 
symmetric naked singularity formation due to the accretion of a null dust fluid. Secondly, the spectra 
of particle creation will provide important informations in order to identity what the products of 
collision are in experiments/observations. However, there is a fundamental problem in estimating the 
spectrum of created particles {i.e., the Bogoliubov coefficients) when the singularity is globally naked: 
we do not know how to impose boundary conditions on a quantum field at the singularity. Thus, 
to avoid such an ambiguity we adopt the model first considered in [12]. This model describes the 
formation of a marginally naked singularity, in which the singularity is observable only within the 
event horizon (so, it violates only the strong version of the Cosmic Censorship Hypothesis [TO])- In 
this spacetime the Cauchy horizon and event horizon coincide, and therefore one has to impose no 
boundary conditions at the singularity. 

This paper is organized as follow. In Sec. [2l we introduce the homothetic or self-similar Vaidya 
solution as the background spacetime describing the formation of the marginally naked singularity. In 
Sec. [3l we consider the creation of massless scalar particles, and find that the spectrum can be recast 
in quasi-thermal form. In Sec. [H we investigate the physical consequences of the result obtained in 
the preceding section, such as the implications of a modified temperature. The final section is devoted 
to conclusion. We use the Planck units unless otherwise noted, in which c = Gd = h = ks = 1- 

think of this point. 
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2 Background spacetime 



2.1 Homothetic Vaidya collapse 

We consider a d-dimensional {d > 4) spherically symmetric spacetime of which hne element is given 
by 

where dO^_2 is the hne element of a unit {d — 2)-sphere. For an arbitrary m(v), the fohowing stress- 
energy tensor solves the Einstein equation, 

Tat = ^^Kh, (2.2) 

where ka := —5^ is an ingoing null vector. Here, we choose the mass function ni{v) as 

' 0, v<0 (I) 

m{v) = < fMv'^-^, 0<v <vo (II) , (2.3) 
^ M, v>vo (III) 

where fx, vq, and M are all positive constants. Then, the metric (|2.ip describes the gravitational 
collapse of a null dust fluid or pressureless radiation flowing toward the center from the asymptotic 
region. The above choice of m{v) results in the homothety or self-similarity of the solution. It is 
known that a central singularity appears at {v, r) = (0, 0), and whether the singularity is naked or not 
depends on the parameters. Usually, the mass parameter M of the Schwarzschild region III is chosen 
as M = ^Vq~^ ^ which results in that the local event horizon in the Vaidya region lj§ is connected 
to that in region III. For now, we do not impose such a condition, and three parameters {fi,VQ,M) 
are left to be independent. We impose a different condition to realize a marginally naked singularity 
in the next subsection. It is noted that the physical mass defined in the asymptotic region is given 

MsH = il^|l^M , (2.4) 

where 0.d~2 = 27r('^~^)/^/r[(d - l)/2] is the volume of the unit {d — 2)-sphere. 

In region II we introduce a new time coordinate (" and similarity coordinate z as 

C ■■= In (-) , z:=-, (2.5) 
\vij r 

where vi is a positive constant^ Then, the line element in region II is written as 

dsl = vfe'^ (^dC^ - |dCdz + ^,dnU^ , (2.6) 

^We mean by the local event horizon the event horizon when the spacetime is entirely filled with the null dust. 
^This parameter vi is introduced just for the dimensions of line element to be correct, but any results in this paper 
do not depend on vi. 
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where 



h{z) := 2fiz'^-'^ - z + 2 . (2.7) 

Furthermore, if we introduce 



the hne element (j2.6p is written in double null form, 

d4 = z;iV^ 0dCd77 + ldJ7i2) ■ (2.9) 

Here, z is understood as an implicit function of ^ and i] via Eq. (|2.8p . By a simple argument (see, 
e.g., [E]), one can show that the singularity appearing at {v,r) = (0,0) is naked when the algebraic 
equation h{z) = [see Fig. [1] (a)] has two distinct positive roots, which we denote by and z+ 
{0 < Z- < z^). For example, in d = 4 the roots can be obtained explicitly as z-t = (1 it ^/l — 16^) /4/j, 
for < ^ < 1/16. In general dimensions, the two positive roots exist if the accretion parameter is 
in the range of 

°<^<^--= [2(d-2)]^-^ ' (2-10) 
and in this case z± is restricted as 

2(d-2) 

2 < z^ < Zc < z+ , Zc := ■ (2-11) 

It turns out that z = Z- and z = z+ corresponds to the Cauchy horizon and event horizon, respectively. 
A conformal diagram of the Vaidya solution with fi satisfying condition ()2.10p is drawn in Fig. [T] (b) . 

2.2 Connecting Cauchy horizon to outer event horizon 

In order to obtain a marginally naked-singular spacetime, the Cauchy horizon in region II, z = z_, is 
connected to the event horizon in the region III, r = rn ■= (2M)^/('^~^) , at the end of accretion of the 
null fluid, V = vq. That is, we impose the following relation among three parameters {fi,VQ, M), 

(2.12) 



(2M)i/("!-3) 



In this case, one can show that mass parameter M is greater than the final quasi-local mass of the 
Vaidya solution, fiv^"^, as follow. From the matching condition (|2.12p and the conditions of z_ < z^ 
and n < He, given in Eqs. (j2.10p and (j2.1ip respectively, we have 



^^-3 = J_ . < =^ < 2Mzf 3 < 2^,zf-3 (2.13) 



2/2 M " M 
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(a) (b) 

Figure 1: (a) Two roots of h{z) = 2fiz'^~'^ — {z — 2) = determines the Cauchy horizon z = z- and event 
horizon z = z^. In the limit of fi ^ fic the two roots degenerate, i.e., z± = z^ ■ (b) A conformal diagram 
of the spacetime entirely filled with the null dust fluid. 

With the help of exphcit expressions of /ic and Zc, Eq. (|2.13|) leads to 

M - > (d - 3)/xi;^-^ > , (0 < V < ^^c) , (2.14) 

which is the desired inequality. Thus, the matching condition ()2.12p and the marginally naked sin- 
gularity is realized, e.g., by infalling a shell of mass M — fJ-VQ~^ toward the center at the moment of 

V = Vq. 

3 Particle creation 

3.1 Near-horizon behavior of massless scalar field 

The Klein-Gordon equation for a massless scalar field is given by 

□(/. = -^da [V^g'^'dbcp) = . (3.1) 

In the rest of this subsection we consider this equation in region IL Separation of variables is possible 
in the (C, 2;)-coordinates (j2.6p by writing 

<P = e<l''Q{z)Yd-2 , (3.2) 

where a is a constant of separation. Yii-2 is a harmonic function on the unit [d — 2)-sphere that 
satisfies 

[Arf_2+£(£ + d-3)]yd_2 = , £ = 0,1,2,..., (3.3) 
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r = 

singularity / + 




Figure 2: A conformal diagram of the spacetime considered in this paper. The Cauchy horizon {z = Z-) 
in the Vaidya region II is connected to the event horizon (r = ru) in the Schwarzschild region \\\ ztv = vq 
by imposing condition ()2.12p . 



where lS.d-2 is the Laplacian on the sphere. The equation for Q{z) takes form of an eigenvalue problem 
for o", 



zhdlQ+ 2^iz'^~^ + {d-A)z + a -{d-A) d,Q + 



Q = 



(3.4) 



where we have partially used the explicit expression of h{z), Eq. ([27 

While it seems difficult to obtain an analytic solution to Eq. (j3.4p globally, it is enough for our 
aim to solve it near the event horizon, z = z_, and the boundary between the regions I and II, z = 0. 
To consider Eq. (13. 4p near z = z_, it is useful to introduce h{z) by 



h{z) = {z- z^)h{z) , h{z^) / 



(3.5) 



where the non- vanishing of h{z^) corresponds to that h has a simple zero at Z-. Then, one can show 
that Eq. (j3.4p has power solutions near z = z_, 

(T + {d-2) 



Q cii {z- z_)' 



, 



■ ^ 2(d- 2) - {d-3)z- ' 
where we have used I'Hopital's rule to eliminate h{z-) in the power, 

Z — Z- Z^Z^ Z- 

On the other hand, Eq. (j3.4p has power solutions near the I-II boundary, z = 0, given by 



Q^z'-^, 



d-2 



a 
2 



(3.6) 



(3.7) 



(3.^ 
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3.2 Quasi-thermal spectrum 

Now, we quantize the massless scalar field of which classical behaviors in region II were investigated 
above, and calculate the spectrum of particle creation largely according to Hiscock et al. [12]. 

A mode function which defines the positive energy particles {dtPuj = —iujpuj, u; > 0) at the future 
null infinity /"*" is written as 

_ N P^jr) 



where is a normalization constant, and u is the retarded time in the region III: 

u:=v-2n, n= j-rjz:^ ■ ^-^0 

J 1 - (rH/ry 

Hereafter, we omit the dependence of mode functions on the spherical coordinates, i.e., 1^-2; since it 
plays no important role in our analysis. Pi,j{t) is a function that approaches the unity at the asymptotic 
region, \\m.j.^ao Puiif) = 1- One the other hand, the mode function which defines the particle at the 
past null infinity is written as 



where N' is a normalization constant and limr^oo F^' (r) = 1 is assumed. 

Since both and fi^' are complete sets, one of them can be expressed as a sum over the frequencies 
of the other. Namely, 



oo 



Pu= I {ctuju,' L' + f^LJLu' , (3.12) 



where a^iui' and fS^^ui' are constants (the Bogoliubov coefficients) and the bar represents the complex 
conjugate. If we assume that the quantum state is in the vacuum |0) defined by f^^r at /~ and observe 
the particles defined by p^j at , we have the number density given by 

roc 

{0\NJ0) = / |/3,^.|'da;' , (3.13) 



where N^^ is the particle number operator corresponding to Pi^. Thus, what to do in order to estimate 
the number density spectrum p.l3p is to propagate to /~ and calculate the Bogoliubov coefficient. 

By tracing p^ along the event horizon to the II-III boundary, v = vq, we find 



W - — - , (3.14) 



where Kq is the product of the normalization constant, a scattering amplitude for the static geometry 
in region III, and the asymptotic limit of the factor P^^/r^*^"^)/^ at r — )• ru- Writing r^, in terms of r, 
Eq. ()3.14|) is written at the leading order, 

p^^^[ 1 , (3.15) 
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where unimportant factors are absorbed in Ki, and k is the surface gravity, 



1 ^ /'rH\d.~3 



K := -dr I 1 — ( — 



2 \ \ r 



d — 3 , , 

(3.16) 



Furthermore, writing r in terms of z, Eq. (j3.15p becomes 

p^^^i^z^-zr'"' , (3.17) 



where unimportant factors are absorbed in K2. Now, we would hke to relate p^j given by Eq. (13.17P 
with the the general solution given by Eq. (|3.6p . Clearly, mode function (j3.17p corresponds to the 
solution with the power in Eq. (|3.6p . Thus, we have a relation between oj and a 

, , . 2(d - 2) - (d - 3)z_ 
a = -{d-2)+ixu; , x ■= Hr^ —■ (3.18) 

From Eq. (|3.8p . the general solution near the I-II boundary, v = 0, can be written as 

„,<T/2+(<i-2)/2 

Ae^l'^z-'^l'' + Be-^/2z('^-2)/2 = ^^-/2 + ^ (3.19) 

where A and i? are constants, and we have set v\ to unity for notational simplicity. Thus, if we trace 
the wave (j3.17p toward the neighborhood of f = 0, and put it in the form of general solution (j3.19p . 
it should take form of 



where we have expressed a in terms of x through Eq. p.lSp . The first term in Eq. ()3.20p corresponds 
to the ingoing wave at z; = 0, while the second term corresponds to the wave scattered to /~ near 
the singularity. Thus, and are the transmission and scattering coefficients, respectively, for the 
dynamical region II. They are related each other by the probability conservation — = 1. 

Tracing the wave (|3.20p to /" via f = and then via the center r = 0, mode function at |u| ~ 
is found to be 



Pi^ — { ^iuX^nv ' (3.21) 



where a phase shift of the transmitted wave is absorbed in A'^. 

Now, we are ready to calculate the Bogoliubov coefficients with using Eq. p.2ip . The coefficients 
are given by reversing Eq. (j3.12p . 

/oo rOD 
e'^'y^dv , = Va7r('^-2)/2 / ^-^-'y^^y . (3.22) 
-00 J —00 
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Substituting in Eq. ()3.2ip into Eq. ()3.22p , we have 



a^^> = K2\ - {A'^J- + B^J+) , P^^, = K2\ - {KJ+ + B^J.) , (3.23) 

where 

J± := / exp (zbiw't! + icjx In df . (3.24) 
Jo 

This integration is nothing but the one evaluated by Hawking [3]. Thus, the resultant expressions of 
(ctujuj' , l^uiui' ) are given as 



1..' (-ALe^'^>^/2 ^ 5^e--x/2^ ^ ^ (^^;^e— x/2 _ i5^e-x/2^ ^ (3.25) 

where 



-juiui' := iK2\ — r(l + iwx) exp [-(1 + iwx) Inw'] . (3.26) 



From Eq. p.25p . we obtain a relation 



2 I — e'^^^ I ^ 



Conservation of probability may be expressed as 

{\a^co'\^ -\f3^^'\^)dio' = T^ , (3.28) 



00 



where is the transmission coefficient for the wave propagating backward from through the static 
geometry in region III and the low-curvature portion of region IL Now, we obtain the spectrum (j3.13p 
with using Eqs. (13:2711 and ([5:28]) . 

(0|^-|0) = ,./tJ_i ' ^^*^^=ln[l + (e2-x-r)/|l-ae-x|2] ' (^.29) 

where we have introduced an frequency-dependent effective temperature to put the spectrum in the 
quasi-thermal form. 



4 Modified temperature and lifetime 

In order to see the physical meaning of the effective temperature p.29p . let us neglect the scattering 
near the singularity {i.e., take the limit of C^^ — t- 0) and define a temperature in this limit, 

T := hm r,ff = = f -Th , Th = ^ , (4.1) 

c^^o 2ttx [d — 3)[Zc — Z-) AirrH 
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0.2 0.4 0.6 0.8 1 

(b) 

Figure 3: (a) Accretion parameter fi v.s. temperature T normalized by the Hawking temperature Th [see 
Eq. (j4.ip ] for d = 4 (solid), 5 (dashed), and 10 (dotted). The temperature diverges in the limit of — Hc- 
(b) Accretion parameter v.s. evaporation time normalized by that determined by the Hawking temperature 
[see Eq. 



where Th = k/2tt is the Hawking temperature. Reviving natural constants (c, h,kB, G^), Th is written 
in term of Mbh, 

where Mp := {h'^~^ /Gdc'^~^y^^^~'^^ is the d-dimensional Planck mass. Numerical values of Hawking 
temperature (j4.2p when Mp(? = 1 TeV and Mbh(? = 5 TeV are given in Table [TJ Since T is an 
increasing function of z_, its minimum value is realized in the limit of — )• 2 or equivalently ;U — )• 0. 
Namely, T > lim2__j.2 T = Th- Thus, the temperature is always higher than the Hawking temperature. 
One can also show that in the limit that the two homothetic horizons coalesce (realized by z± — ?• Zc or 
jj, He) the temperature diverges, lim^ -j.^^ T = +oo. The /i-dependence of T is shown in Fig. [3] (a). 

If we take into account the backreactions of particle creation, the modification of temperature 
affects the evolution of the spacetime, e.g., the lifetime or evaporation time. The mass decreasing 



Table 1: Hawking temperature Th and its corresponding evaporation time th when Mpc^ = 1 TeV 
and Mbhc^ = 5 TeV. For reference, Schwarzschild radius rn is given in the unit of fermi (1 fm := 
10"^^ m). 



d 


5 


6 


7 


8 


9 


10 


11 


UbTh (GeV) 


77.2 


179 


282 


380 


470 


553 


629 


TH (10-26 sec) 


220 


24.5 


5.31 


1.44 


0.418 


0.123 


0.0357 


rn (10-4 fm) 


4.07 


2.64 


2.23 


2.07 


2.01 


1.99 


2.00 
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rate is proportional to the luminosity, which is the product of horizon area Ah and flux oc T'^ (i.e., 
Stephen-Boltzmann law in d-dimensional spacetimes) , 

-^ocyi^rVM-^-^), (4.3) 

where we have used Ah fx r^~^ cx M^j'^^^-'^'-'^"^'' and T oc M^^^'^'~^\ Integrating (|4.3p . the evaporation 
time T is estimated as 

r oc oc r-('^-i) . (4.4) 

Taking into account that Eq. (|4.ip holds in our case, we obtain 

' ^'-'^^;'-'"^ ]"\h, (4.5) 
where th is the evaporation time of a black hole of temperature Th, 

~ Mpc"^ V Mp ) [ (d - 2)2 ) 



(2(i-l)/2 



(d - 3)'^(d - i)2r[(d + 2)/2]C(d) ■ 

(4.6) 

Here, ^ is the Riemann zeta-function. The numerical values of th are given in Table [ij^ Thus, the 
lifetime is always shorter than that of an ordinary black hole of the same mass, and even can be 
arbitrarily short in the limit that the two homothetic horizons coalesce. The /z-dependence of r is 
shown in Fig. [3] (b). 



5 Conclusion 

Motivated by the possible production of naked singularities by high-energy particle collisions in the 
context of TeV-scale gravity, we have investigated the particle creation by naked singularities in 
higher dimensions. For simplicity, the background spacetime has been modeled by the d-dimensional 
spherically symmetric collapsing solution, given by Eqs. (|2.1|) and (|2.3|) . In addition, the Cauchy 
horizon in the Vaidya region is matched to the event horizon in the outer region by imposing a 
particular condition ()2.12p . This matching enables us to calculate the Bogoliubov coefficients without 
the ambiguity of boundary conditions at the singularity. We have found that the spectrum of the 
created particles deviates from that of the usual Hawking radiation due to the scattering near the 
singularity. However, the spectrum can be recast in the quasi-thermal form with a frequency-dependent 

''Here, th is the evaporation time of a black hole in a minimal scenario. That is, we assume that only gravitons, of 
which polarization degrees of freedom is (d — 3)d/2 in d-dimcnsions, propagate in the bulk (see the first paper in Ref. 
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temperature, Eq. (|3.29p . Neglecting the scattering near the singularity, the temperature is defined as 
Eq. ()4.ip and compared with the Hawking temperature. The temperature is always higher than the 
Hawking temperature of a same-mass black hole, and therefore the evaporation time is always shorter 
than that of the same-mass black hole, in the whole range of parameters (i.e., spacetime dimensions 
d and accretion parameter /i). This is a concrete prediction about the particle creation during naked- 
singularity formation in higher dimensions, made by this paper for the first time (as far as the authors 
know) . 

One might wonder that matching condition (|2.12p choses a particular class of solutions out of more 
general solutions. While that condition is essential to avoid the ambiguity of boundary conditions at 
the singularity, there seems no physical reason why such a class of solutions with the marginally naked 
singularity are preferred. Therefore, a more general case including globally naked singularities should 
be analyzed. The present authors are working in this direction and will report it elsewhere soon jl9j . 

Here, we discuss the relation between the analysis in this paper with that in [TJ]. Paper [ll] 
argued that (in the usual 4-dimensional gravity) if unknown quantum gravitational effects dominate 
near the birth of naked singularity to suppress the particle creation, the net energy radiated amounts 
only to a few Planck energy Mp(4) ~ 10^^ TeV. If this is the case, the particle creation by the naked 
singularity forming in the gravitational collapse of a massive object, e.g., an astrophysical star of mass 
Mq ~ 10^^ TeV, will be a negligible effect. However, we stress that in the TeV-scale gravity scenarios, 
even if the quantum gravitational cutoff works, the particle creation by the naked singularity can be 
a crucial effect. Namely, even if the net energy radiated by the particle creation amounts only to a 
few Planck energy ~ O(TeV) due to the cutoff, the naked singularity produced by particle collisions 
{e.g., in nowadays collider experiments) itself would have the same order of mass at most. Therefore, 
the dynamics of background geometry is significantly affected by the particle creatior^. 

The production of naked singularities (or visible borders of spacetime) in the TeV-scale gravity 
is an intriguing possibility. However, we have not gotten any evidence of the production of a naked 
singularity nor a black hole in particle collisions so far. In recent years, the Numerical General 
Relativity is expanding to higher dimensions. In fact, several results on head-on collisions of black 
holes in higher dimensions have been reported [20]. (See [21] and [22] for the high-energy collisions 
of 4-dimensional black holes with zero and non-zero impact parameters, respectively.) While the 
collision of higher-dimensional black holes would be different from that of Standard Model particles, 
which are confined to our 3-brane, it would be interesting to look for any signs or symptoms of 
naked-singularity /visible-border formation in numerical experiments. 

""Incidentally, we should repeat here that the present model is quite simplified, and therefore the validity of the 
approximations used in this paper should be verified carefully when we discuss the actual process of visible-border 
production. For instance, the visible-border production is expected to be a highly dynamical process, of which typical 
timescale may be given by the light-crossing time through the collision region, and therefore the validity of the geometric- 
optics method and/or the QFT in curved background itself could be questionable. The authors thank an anonymous 
referee for focusing our attention to this point. 
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